Journal: Frontiers in Network Physiology
Article Title: A guide to Whittle maximum likelihood estimator in MATLAB
doi: 10.3389/fnetp.2023.1204757
Figure Lengend Snippet: Whole set of α ^ estimated. On the x -axis, α is the true value of the exponents, the values computed in the generator. On the y -axis, the estimated values of α ^ computed by the corresponding analysis method are located. Red curves represent true alpha values α , and blue curves represent the estimated alpha values α ^ . The first column (A,C,E) corresponds to the signals generated via the Cholesky method, and the second column (B,D,F) corresponds to the signals generated via ARFIMA filtering. The first row (A,B) presents the α ^ values computed using fGn-based Whittle’s maximum likelihood estimator, the second row (C,D) presents those computed using ARFIMA-based Whittle’s maximum likelihood estimator, and the third row (E,F) presents the α ^ values computed using DFA.
Article Snippet: .................................................................. MATLAB code 6 : Whittle’s log-likelihood function with fGn theoretical PSD function lwHfgn = WLLFfgn (H,w,P,N) Tp = sin (pi*H)*gamma ((2*H)+1)*(abs(w).^(1-(2*H))); c = sum(P)/sum (Tp); T = c *Tp; lwHfgn=(2/N)*sum (log(T)+(P./T)); .................................................................. MATLAB code 7 : Whittle’s log-likelihood function with ARFIMA (0,d,0) theoretical PSD function lwHarf = WLLFarf (H,w,P,N) d = H-0.5; Tp=(1/(2*pi))*(2*sin (w/2)).^-(2*d); c = sum(P)/sum (Tp); T = c *Tp; lwHarf=(2/N)*sum (log(T)+(P./T)); .................................................................. In the first line, the function is used to declare the functions WLLFfgn and WLLFarf .
Techniques: Generated